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A  common  way  to  model  lithium-ion  batteries  is  to  apply  equivalent  circuit  (EC)  models.  In  this  work 
two  different  EC  models  are  build  up  and  parameterized  for  a  commercial  6.5  Ah  high-power  lithium-ion 
cell.  Measured  impedance  spectroscopy  data  depending  on  temperature  and  state  of  charge  (SOC)  are 
used  for  parameter  estimation. 

The  first  EC  model  consists  of  an  ohmic  resistor  (R),  an  inductor  (L)  and  three  RC-elements  (a  parallel 
connection  of  a  capacitor  (C)  and  a  resistor).  The  second  EC  model  consists  of  one  R,  one  L,  two  Zarc  ele¬ 
ments  and  a  Warburg  element.  The  estimated  parameters  were  used  to  develop  two  empirical  electrical 
cell  models  which  are  able  to  predict  the  voltage  of  the  cells  depending  on  current,  temperature  and  SOC. 
Hereby  the  internal  cell  resistance  Ri  is  based  on  the  EC  models  and  a  Butler-Volmer  adjustment.  Both 
approaches  were  validated  by  current  profiles,  which  cover  typical  automotive  applications  to  prove  the 
model  performance  at  low  temperatures  and  high  dynamic  operation.  An  accurate  voltage  prediction 
could  be  realized  with  both  EC  models.  The  second,  more  complex,  model  is  able  to  predict  cell  voltage 
more  precisely,  but  at  the  expense  of  up  to  four  times  higher  computational  effort. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  order  to  predict  the  performance  of  a  battery  cell  depending 
on  operation  conditions  and  load,  simulations  based  on  an  equiv¬ 
alent  circuit  (EC)  model  are  widely  used  [1,2].  In  this  work  two  EC 
models  with  different  complexities  are  presented  and  compared 
in  terms  of  simulation  performance  and  effort.  The  generalized 
method  to  identify  these  models  is  presented  in  Fig.  1. 

In  the  first  step  the  parameters  of  the  model  have  to  be  esti¬ 
mated.  Therefore  impedance  spectra  gained  by  the  EC  are  fitted 
to  experimental  measured  data  through  a  numerical  optimization. 
Afterwards  the  estimated  parameters  are  implemented  into  a  cell 
model.  Finally,  the  model  is  validated  by  comparison  to  measured 
data. 

A  precise  prediction  of  voltage  allows  testing  battery  cells 
under  different  conditions  just  by  simulation  and  therefore  reduces 
experimental  efforts  and  time.  Moreover  simulations  allow  an  opti¬ 
mization  of  the  battery  and  in  the  next  step  of  drive  trains  for 
upcoming  hybrid  and  electrical  vehicles. 
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In  a  previous  work  [3]  impedance  spectra  were  recorded  for  a 
high-power  6.5  Ah  lithium-ion  cell  in  a  temperature  range  from 
-30  °C  to  50  °C  for  the  entire  state  of  charge  (SOC)  range.  These 
spectra  contain  all  information  about  cell  impedance  depend¬ 
ing  on  operating  conditions  and  allow  therefore  a  prediction  of 
the  cell  performance  under  different  loads  and  scenarios.  Based 
on  these  spectra  both  EC  models  are  set,  verified,  validated  and 
compared  in  terms  of  computation  time  and  accuracy  of  voltage 
prediction. 

2.  Theory 

2.1.  Equivalent  circuit  with  RC-elements 

The  first  applied  EC  to  model  cell  voltage  was  built  up  in  a  most 
simple  way  in  order  to  have  a  basic  model.  Aim  of  this  approach 
is  to  gain  a  very  fast  computing  model  and  to  get  a  benchmark  for 
an  advanced  EC  model.  The  EC  consists  of  one  element  for  every 
section  of  the  ideal  impedance  spectrum  of  lithium-ion  cells  and  is 
displayed  in  Fig.  2. 

The  impedance  behaviour  of  the  cell  at  very  high  frequencies 
is  represented  by  an  ideal  inductor  L  and  the  ohmic  resistance  by 
a  resistor  R.  The  following  electrochemical  processes  at  the  elec¬ 
trodes  can  be  represented  in  a  most  simple  way  by  RC-elements, 
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Fig.  1.  Approach  of  system  identification  for  an  EC  model. 


a  parallel  connection  of  a  capacitor  C  and  a  resistor  R.  Thus  the 
impedance  of  the  system  can  be  described  as  the  sum  of  its  ele¬ 
ments: 

-RC  =JOjL  +  R*  +  1  Jr~2  +jojC2  +  1  Jr~3  +JojC3  +  1/R4  +JcoC~4  ' 

This  approach  is  a  state  of  the  art  EC  [4-7]  and  has  in  total  eight 
parameters  which  have  to  be  estimated  and  fitted  by  an  optimiza¬ 
tion  algorithm. 

2.2.  Equivalent  circuit  with  Zarc  and  Warburg  element 

The  goal  of  the  second  EC  is  a  most  accurately  and  meaningful 
reproduction  of  the  measured  impedance  data.  An  approach  with 
Zarc  and  Warburg  elements  suggested  by  several  authors  [8-10]  is 
shown  in  Fig.  3. 

2.2.1.  The  Zarc  element 

Zarc  elements  are  a  parallel  connection  of  a  resistor  R  and  a 
constant-phase-element  (CPE)  [11].  CPEs  consist  of  a  generalized 


Fig.  3.  Ideal  impedance  spectrum  of  a  lithium-ion  cell  and  an  EC  with  Zarc  and 
Warburg  element. 


capacity  0  and  depression  factor  and  are  defined  by 

ZcpE'S5v 

A  parallel  connection  with  a  resistor  leads  to  Eq.  (3): 

z  i 

-Zarc  VR+(ja>/e' 


(2) 

(3) 


Comparing  Eq.  (3)  with  the  RC-element  of  Eq.  (1)  both  equations 
are  equal,  if  is  left  aside. 

This  depression  factor  iff,  valid  between  zero  and  one,  is  respon¬ 
sible  for  the  depression  of  the  semi-circle  in  a  Nyquist  plot  as  Fig.  4 
illustrates.  For  \//  =  0  the  Zarc  element  represents  only  an  ohmic 
resistance  and  for  \fr  =  1  the  ideal  semi-circle  of  a  RC-element.  The 
advantage  of  Zarc  elements  is  their  potential  to  reproduce  the 
depression  of  semi-circles  in  impedance  spectra,  which  are  char¬ 
acteristic  for  lithium-ion  cells.  However,  the  disadvantage  of  Zarc 
elements  is  a  missing  Laplace  transformation.  Thus  it  is  not  possible 
to  transfer  the  parameter  sets  from  the  frequency  domain  into  the 
time  domain  without  approximations  [12]. 

Buller  [13]  introduced  in  his  thesis  a  way  to  approximate  Zarc 
elements  by  a  variable  number  of  RC-elements  as  illustrated  in 
Fig.  5. 

The  principle  of  this  approximation  is  based  on  an  odd  num¬ 
ber  of  semi-circles  whose  radii  are  increasing  towards  the  medium 
semi-circle.  The  resistance  and  capacity  of  each  RC-element  can  be 
computed  by  optimization  factors  which  are  listed  in  the  literature 
[  1 3  ].  In  Fig.  6,  the  result  of  an  approximation  with  five  RC-elements, 


x  1CT4 


ip  =  1  - ip  =  0.8  - tp  =  0.6  — ■ —  ip  =  0.4 


Fig.  2.  Ideal  impedance  spectrum  of  a  lithium-ion  cell  and  an  EC  with  three  RC-  Fig.  4.  Impact  of  depression  factor  f  on  the  Nyquist  plot  of  the  Zarc  element 
elements.  {R  =  l  m£2,  C=  1  kF). 
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Fig.  5.  Approximation  of  the  Zarc  element  with  five  RC-elements. 
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Fig.  6.  Nyquist  plot  of  a  Zarc  approximation  with  five  RC-elements,  the  dash-dotted  Fig.  7.  Nyquist  plot  of  Warburg  element,  the  solid  line  represents  the  ideal  Warburg 

line  represents  the  ideal  Zarc  element,  the  dashed  line  represents  the  circuit  with  element  and  the  dashed  line  represents  the  R-RC  circuit, 

five  RC-elements  and  the  solid  line  represents  the  circuit  with  just  one  RC-element 
(R  =  10  mQ,  C= 2 1<F  and  f  =  0.9). 


as  used  in  this  work,  in  comparison  to  an  approximation  with  just 
one  RC-element  is  shown. 

Especially  for  high  values  of  ^  the  approximation  matches  the 
ideal  curve  very  accurate  and  is  therefore  valid. 


2.2.2.  The  Warburg  element 

Analogous  to  the  Zarc  element,  an  approximation  is  required  for 
the  Warburg  element.  A  Warburg  element  is  defined  by 


7 

— w 


Rtanh (yjw0) 

VJojd 


(4) 


and  shows  a  linear  slope  with  an  angle  of  45°  at  high  frequencies. 
This  linear  slope  predestinates  Warburg  elements  to  represent  the 
diffusion  process  of  Fig.  3. 

A  simple  way  to  approximate  the  Warburg  element  is  a  series 
connection  of  a  resistor  and  a  RC-element  as  illustrated  in  Fig.  7. 
An  approach  applied  by  Handschuh  [14].  The  impedance  for  this 
connection  is 


—W  —  Ri  + 


1 

1/^2  +  jti>C2 


(5) 


For  the  estimation  of  the  parameters  R 1 ,  R2  and  C2,  the  frequency  coq 
where  the  spectrum  reaches  its  minimum  value  is  required.  Hand¬ 
schuh  [14]  showed  that  the  value  of  co0  is  independent  of  R  and  can 
be  approximated  well  by 


2.53 

w°  =  ^r- 


(6) 


Furthermore  it  was  shown  that  the  minimum  value  of  the  spec¬ 
tra  depends  only  on  the  resistance 


By  using  this  equation,  the  parameters  can  be  determined: 


«.  =  §«• 

(8) 

*  =  §«, 

O) 

1 

6  1 

c2  =  „ 

CO0R2 

5  coqR' 

(10) 

In  order  to  get  an  approximated  slope  of  45°,  R 1  has  to  be  separated 
in  additional  two  RC-elements  with  following  parameters: 


«...  -  n«- 

(ii) 

(12) 

Ci.1  =  ^qC2; 

(13) 

Cl, 2  =  gC2. 

(14) 

Finally,  this  approximation  allows  to  representing  the  slope  with 
high  accuracy  as  demonstrated  in  Fig.  8. 

Since  for  this  work  just  the  linear  part  of  the  Warburg  element 
is  required,  the  last  approximation  is  done  for  the  intersection  of  a 
tangent  which  runs  through  the  point  of  origin  and  the  RC-element 
by  their  geometrical  relation: 


c3- 


2.42  (1  -  a/05) 


c2. 


(15) 

(16) 


Z"  ■ 

min 


By  using  these  approximations  it  is  possible  to  represent  the  War¬ 
burg  element  accurately.  Thus  the  impedance  of  the  advanced  EC 
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Fig.  8.  Nyquist  plot  of  a  Warburg  approximation  with  three  RC-elements,  the  red 
line  represents  the  ideal  Warburg  element  and  the  blue  line  represents  the  3-RC 
circuit.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  the  article.) 

can  be  computed  by  the  equation  corresponding  to  Fig.  3: 

1  1 

Zzw-J<0L  +  Rl  +  —  -  ^  ^ 

+  Rtanh (Vjw&i) 

V^4 


was  used  for  every  spectrum.  For  in  total  119  spectra,  the  parame¬ 
ters  were  estimated  and  interpolated  to  a  characteristic  map.  Since 
for  low  temperatures  no  diffusion  section  was  recorded  by  the 
electrochemical  impedance  spectroscopy  (EIS)  because  of  too  low 
frequencies,  an  additional  algorithm  was  required.  This  algorithm 
was  able  to  detect  the  existence  of  a  diffusion  section  in  the  mea¬ 
sured  data  and  was  setting  R3  and  C3  to  zero  if  no  diffusion  process 
was  recorded. 

The  method  of  optimization  is  illustrated  in  Fig.  9  by  plotting 
measured  and  fitted  curves  for  an  exemplary  SOC  of  60%  at  different 
temperatures  in  Nyquist  plots.  The  general  behaviour  is  similar  for 
all  temperatures.  The  simple  EC  is  able  to  approximate  the  general 
behaviour  of  the  curves,  but  fails  to  reproduce  the  spectra  accu¬ 
rately.  Problematic  are  the  loops  of  the  measured  spectra  which 
are  depressed  semi-circles  and  therefore  difficult  to  fit  by  ideal  ones 
and  the  almost  linear  diffusion  part.  Therefore  large  deviations  can 
be  seen  especially  in  the  middle  of  the  first  and  second  RC-element 
and  for  the  linear  slope  at  low  frequencies.  Since  very  high  frequen¬ 
cies  have  no  practical  relevance  for  the  simulation  of  batteries,  the 
inductivity  was  neglected  and  is  not  represented. 

3.2.  Equivalent  circuit  with  Zarc  and  Warburg  element 


In  total  this  EC  has  ten  parameters,  two  parameters  more  than  the 
first  EC,  increasing  the  computation  time  but  also  the  accuracy  of 
the  simulation. 

3.  Parameter  estimation 

3.1.  Equivalent  circuit  with  RC-elements 

In  this  section,  the  parameters  of  the  introduced  ECs  are  esti¬ 
mated  by  numerical  optimization.  For  this  purpose  a  nonlinear 
curve-fitting  optimization  based  on  the  method  of  least-squares 


The  same  optimization  approach  was  applied  to  the  second  EC 
with  Zarc  and  Warburg  elements.  The  EC  with  estimated  parame¬ 
ters  was  again  compared  to  measured  data  and  the  result  is  given 
in  Fig.  10. 

It  is  obvious  that  this  EC  represents  the  measured  data  with  a 
higher  accuracy.  Both,  the  semi-circles  as  well  as  the  diffusion  part 
are  approximated  very  precisely  by  the  Zarc  and  Warburg  elements. 
Especially,  the  depressed  semi-circles  are  represented  very  well 
by  the  depression  factor  t/c  Flowever,  small  variations  can  still  be 
seen  for  elevated  temperatures  at  the  change  from  the  second  Zarc 
element  to  the  Warburg  element. 


Temperatur  =  -30°C 


x  i  o-3  T emperatur  =  4°C 


Z’/Q  xIO"3 


Temperatur  =  -24°C 


x  1 0-4  T  emperatur  =  24°C 


Temperatur  =  50°C 


Fig.  9.  Comparison  of  measured  and  simulated  EIS  data  by  the  EC  with  RC-elements  for  different  temperatures  at  60%  SOC,  the  solid  line  represents  the  measured  spectra 
and  the  dashed  line  represents  the  simulated  spectra. 
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Fig.  10.  Comparison  of  measured  and  simulated  EIS  data  by  the  EC  with  Zarc  and  Warburg  elements  for  different  temperatures  at  60%  SOC,  the  solid  line  represents  the 
measured  spectra  and  the  dashed  line  represents  the  simulated  spectra. 


4.  Modelling 


4.2.  Equivalent  circuit  with  RC-elements 


The  estimated  parameters  were  used  to  build  cell  models.  Input 
parameters  of  this  model  were  current  I,  cell  capacity  C,  cell  tem¬ 
perature  T  and  initial  state  of  charge  SOQ=0.  By  integrating  the 
cell  current,  the  actual  SOC  could  be  determined  continuously.  Fur¬ 
thermore  the  open-circuit  voltage  (OCV)  curve  of  the  cell,  which 
was  shown  in  [1],  was  approximated  as  a  function  of  SOC  by  a 
lookup  table.  By  adding  the  overvoltages  q  of  the  system,  which 
are  provided  in  form  of  the  EC  parameters,  the  actual  cell  voltage 
was  computed  and  visualized. 

4.1.  Butler-Volmer  adjustment 

Since  no  DC-offset  was  used  for  the  EIS  measurements  the 
effect  of  current  on  cell  overvoltage  was  unknown.  However,  this 
effect  is  important  for  the  charge  transfer  and  is  defined  by  the 
Butler-Volmer  (BV)  equation: 

i  =  iQ  .  ( eazFIRT(E-Eeq )  _  e-((l-a)zF/Kr)(£-£eq))  (18) 

where  i  is  the  electrode  current  density,  i0  is  the  exchange  cur¬ 
rent  density,  a  is  the  symmetry  factor,  z  is  the  number  of  electrons, 
E  is  the  electrode  potential,  Eeq  is  the  equilibrium  potential,  F  is 
the  Faraday  constant,  R  is  the  universal  gas  constant  and  T  is  the 
temperature. 

In  order  to  consider  the  impact  of  current,  a  comparison  was 
done  between  simulated  (EC  with  RC-elements)  and  measured 
overvoltage  for  18  s  long  pulses  at  60%  SOC  at  different  temper¬ 
atures.  The  differences  between  both  curves  were  used  to  generate 
a  lookup  table  with  adjustment  factors  for  the  overvoltage  as  a 
function  of  temperature  and  current.  Moreover  it  was  observed 
that  positive  and  negative  current  pulses  cause  very  similar  voltage 
responses.  In  consequence  the  symmetry  factor  a  was  set  to  0.5. 


Three  parts  are  necessary  for  modelling  the  actual  cell  voltage. 
The  first  one  is  a  block  for  the  calculation  of  the  actual  OCV.  By 
integrating  the  current,  the  actual  SOC  is  computed  and  recalcu¬ 
lated  in  the  open-circuit  voltage  by  a  lookup  table.  The  second  one, 
the  ohmic  resistance  depending  on  the  actual  /,  T  and  SOC  is  approx¬ 
imated  by  a  lookup  table.  The  last  part  is  the  overvoltages  by  the 
RC-elements. 

For  implementing  the  estimated  parameters  of  the  RC-element 
in  the  EC  model,  a  correlation  between  input  signal  current  I  and 
output  signal  voltage  U  is  required  for  every  RC-element.  The  cur¬ 
rent  of  this  parallel  connection  is: 

I  =  Ir  +  Ic •  (19) 


The  current  flowing  through  the  resistor  and  capacitor  are  defined 
by 

Ir=j,  (20) 


,  r  3  U 

Ic  =  cw 


(21) 


Including  Eqs.  (20)  and  (21)  in  Eq.  (19),  the  voltage  of  the  RC- 
element  can  be  calculated  by 

In  this  equation  R  and  C  are  functions  of  SOC,  T  and  are 
gained  again  by  lookup  tables.  The  current  I  is  the  time  dependent 
input  signal  and  the  actual  voltage  can  be  calculated  by  numeri¬ 
cal  integration  using  Eq.  (22).  As  introduced  in  Section  4.1  these 
impedances  are  corrected  by  determined  BV  parameters.  Finally, 
the  actual  cell  voltage  is  calculated  by  adding  up  all  three  voltages. 
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x  105 
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Fig.  11.  Stepwise  discharge,  on  top  EC  model  with  RC-elements,  on  bottom  EC  model  with  Zarc  and  Warburg  elements. 


4.3.  Equivalent  circuit  with  Zarc  and  Warburg  element 

The  second  modelling  was  done  similarly  to  the  described  one. 
In  addition,  submodels  based  on  the  equations  of  Section  2.2  for 
the  parameter  estimation  of  the  Zarc  and  Warburg  approximations 
were  built.  Analogous  to  Fig.  1 1  the  submodels  were  connected  to 
a  cell  model  and  the  actual  cell  voltage  could  be  computed.  Due 
to  increased  clarity  in  the  further  text,  the  cell  model  with  RC- 
elements  is  abbreviated  to  RC  model  and  the  cell  model  with  Zarc 
and  Warburg  element  is  abbreviated  to  ZW  model. 

5.  Validation 

Both  models  were  validated  against  different  current  profiles  in 
order  to  check  their  performance  in  a  wide  SOC  and  temperature 
range. 

The  first  validation  was  carried  out  to  test  the  models  in  the 
whole  SOC  range  by  comparing  the  results  of  cell  characterization 
tests  for  simulation  and  measurement.  Moreover,  the  models  were 
tested  and  analysed  for  several  cycling  profiles  and  finally  validated 
against  data  recorded  during  a  hybrid  electrical  vehicle  (HEV)  drive. 

Since  the  computation  time  plays  an  important  role  regarding 
applicability  decisions,  a  summary  of  the  computation  times  for  all 
profiles  is  given  in  Table  1.  The  more  complex  ZW  model  with  in 
total  1 3  RC-elements  is  by  a  factor  up  to  four  slower  than  the  simple 
model  with  just  three  RC-elements. 


Table  1 

Comparison  of  computation  times  for  different  current  profiles. 


Profile 

RC  model  (s) 

ZW  model  (s) 

Factor 

Capacity  test  (~1 0,000  s) 

97 

179 

1.8 

Stepwise  discharge  (~12,000s) 

101 

182 

1.8 

Test  cycle  T  =  -20  °C  (~6100  s) 

40 

96 

2.4 

Test  cycle  T=  12  °C  (~6100s) 

79 

166 

2.1 

Test  cycle  T=50°C  (~6100s) 

58 

231 

4.0 

HEV  cycle  (~2250s) 

41 

105 

2.6 

However,  for  most  profiles  the  computation  time  is  just  doubled 
and  even  for  the  more  complex  model  the  simulation  time  is  always 
below  5%  of  the  real  time. 

5.1.  Stepwise  discharge 

The  first  validations  were  done  for  typical  cell  characterization 
tests  such  as  cell  charging  and  measurement  of  cell  capacity.  The 
result  of  the  cell  discharge  in  steps  of  10%  ASOC  for  measurement 
and  simulation  are  shown  in  Fig.  12.  It  can  be  observed  that  both 
models  were  able  to  predict  cell  voltage  with  a  high  quality  until  a 
SOC  of  1 0%.  From  1 0%  to  0%  the  cell  voltage  is  dropping  down,  since 
the  voltage  reaches  the  end  of  discharge  voltage  limit.  Both  models 
are  not  able  to  simulate  this  phenomenon,  since  the  Butler-Volmer 
coefficient  a  is  probably  varying  at  low  SOC.  In  order  to  compare 
the  results  of  both  models,  the  mean  value  of  the  absolute  error  \s~l 
the  maximum  error  £max  as  well  as  the  sum  of  squared  errors  (SSE) 
between  measured  and  simulated  cell  voltage  are  given  in  addition 
in  all  figures. 

Both  models  show  a  similar  distribution  of  voltage  deviations 
and  have  a  very  comparable  performance  with  a  mean  error  of 
about  8  mV.  No  difference  can  be  seen  for  the  SSE,  indicating  the 
same  quality  of  simulation  for  both  models.  Moreover  a  maximum 
error  of  27%  exists,  but  is  occurring  just  for  a  negligible  number  of 
measurement  points  at  very  low  SOC. 

5.2.  Test  cycle 

The  next  validations  were  done  against  drive  cycles  at  five  dif¬ 
ferent  temperatures  measured  in  the  laboratory.  Only  the  results 
for  12  °C  are  shown  which  are  representative  for  the  observations 
at  other  temperatures.  The  test  cycle  consists  of  a  mix  of  pulses, 
crankings  and  drive  data  for  a  period  of  1.7  h.  Fig.  13  shows  the 
different  grade  of  resolution  for  the  comparison  of  simulation  and 
measurement.  Therefore  in  a  second  figure  (Fig.  14)  just  a  time  slot 
of  the  first  20  min  is  shown. 
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Fig.  12.  Test  cycle  measurement  and  simulation. 


Both  models  showed  again  a  very  satisfying  agreement  with 
the  measured  values.  Especially,  the  zoom  in  Fig.  13  shows  that 
even  the  simple  model  is  able  to  represent  fast  voltage  variations 
within  few  seconds  very  well.  The  error  diagram  shows  very  simi¬ 
lar  values  and  just  a  small  improvement  of  SSE  can  be  seen  for  the 
ZW  model. 

5.3.  Drive  cycle 

The  last  profile  was  measured  during  test  drives  in  South  Africa 
and  is  very  suited  to  test  the  models  in  consideration  of  high  SOC 
variations  (ASOC>30%)  and  temperature  variations. 


By  comparing  the  error  distributions  in  Fig.  14  it  is  obvious  that 
the  RC  model  has  a  wider  variation  and  almost  doubled  SSE  com¬ 
pared  to  the  ZW  model.  Moreover  both  models  underestimate  the 
cell  voltage  most  of  the  time,  especially  at  SOCs  below  60%  which 
is  probably  a  consequence  of  the  BV  adjustment  which  is  only  valid 
for  a  SOC  of  60%.  Even  if  the  voltage  peaks  are  not  simulated  accu¬ 
rately,  both  models  show  in  general  a  very  good  agreement  with  a 
mean  error  below  20  mV  and  predict  therefore  the  cell  performance 
very  accurately. 

A  summary  of  the  simulation  quality  represented  by  the  mean 
error  for  all  tests  is  given  in  Table  2.  Except  one  case  the  more  com¬ 
plex  ZW  model  predicts  the  cell  voltage  with  a  higher  accuracy.  In 
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Fig.  13.  Test  cycle  at  T=  12  °C,  on  top  EC  model  with  RC-elements,  on  bottom  EC  model  with  Zarc  and  Warburg  elements. 
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Fig.  14.  Drive  cycle,  on  top  EC  model  with  RC-elements,  on  bottom  EC  model  with  Zarc  and  Warburg  elements. 


Table  2 

Comparison  of  voltage  deviations  for  different  current  profiles. 


Profile 

RC  model  (mV) 

ZW  model  (mV) 

As  (%) 

Capacity  test 

11.1 

10.5 

5.4 

Stepwise  discharge 

8.3 

7.5 

9.6 

Stepwise  charge 

7.7 

7.2 

6.5 

Test  cycle  T=  -20  °C 

10.6 

10.6 

0 

Test  cycle  T=-8°C 

5.9 

7.7 

-30.5 

Test  cycle  T=  12  °C 

10.3 

9.5 

8.0 

Test  cycle  T=  30  °C 

8.3 

6.7 

19.3 

Test  cycle  T=  50  °C 

7.8 

7.1 

2.6 

HEV  cycle 

18.7 

15.5 

17.1 

general  the  improvement  of  simulation  is  in  the  range  of  5-20%, 
but  is  always  below  3  mV. 

6.  Results  and  discussion 

In  the  last  section,  two  different  EC  models  were  introduced. 
First,  a  simple  model  based  on  an  EC  of  three  RC-elements  and 
second  a  more  complex  model  based  on  Zarc  and  Warburg  ele¬ 
ments.  The  parameters  of  both  ECs  were  estimated  by  a  fitting 
for  a  large  number  of  measured  impedance  spectra.  It  was  shown 
that  the  simple  EC  was  not  able  to  reproduce  the  measured  data 
even  with  optimized  parameters  sufficiently.  In  opposite,  the  more 
advanced  EC  was  reproducing  the  spectra  very  accurate.  Since  for 
Zarc  as  well  as  Warburg  elements  no  Laplace  transformation  exists, 
approximations  were  derived.  These  gained  ECs  were  implemented 
in  cell  models.  Afterwards  both  models  were  validated  against  dif¬ 
ferent  experimental  data  as  cell  characterization  tests,  test  cycles 
and  drive  cycle. 

All  validations  showed  that  both  models  predict  cell  voltage  very 
precisely.  Although  the  simple  EC  was  not  able  to  reproduce  the 
impedance  curves,  voltage  was  predicted  with  a  high  quality.  Com¬ 
parisons  of  both  models  evidenced  an  up  to  20%  higher  mean  error 
for  the  simple  model  by  a  minimum  doubled  computation  speed. 
However,  also  the  mean  error  for  this  model  is  just  in  the  range  of 


a  few  mV.  Therefore  it  can  be  concluded  that  the  impedance  fitting 
plays  a  role  for  the  simulation  quality,  but  the  influence  is  just  small 
and  voltage  can  also  be  predicted  accurately  by  simple  EC  models. 

To  sum  up,  both  EC  models  allow  a  very  precisely  prediction  of 
battery  voltage  and  battery  performance.  Thus  it  is  possible  to  use 
these  battery  models  for  simulations  of  e.g.  consumption,  vehicle 
range,  cooling  requirement  and  battery  lifetime  prognosis. 
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